3.3.71 \(\int \frac {\sec ^{\frac {5}{2}}(a+b \log (c x^n))}{x} \, dx\) [271]

3.3.71.1 Optimal result
3.3.71.2 Mathematica [A] (verified)
3.3.71.3 Rubi [A] (verified)
3.3.71.4 Maple [B] (verified)
3.3.71.5 Fricas [C] (verification not implemented)
3.3.71.6 Sympy [F(-1)]
3.3.71.7 Maxima [F]
3.3.71.8 Giac [F(-1)]
3.3.71.9 Mupad [F(-1)]

3.3.71.1 Optimal result

Integrand size = 19, antiderivative size = 93 \[ \int \frac {\sec ^{\frac {5}{2}}\left (a+b \log \left (c x^n\right )\right )}{x} \, dx=\frac {2 \sqrt {\cos \left (a+b \log \left (c x^n\right )\right )} \operatorname {EllipticF}\left (\frac {1}{2} \left (a+b \log \left (c x^n\right )\right ),2\right ) \sqrt {\sec \left (a+b \log \left (c x^n\right )\right )}}{3 b n}+\frac {2 \sec ^{\frac {3}{2}}\left (a+b \log \left (c x^n\right )\right ) \sin \left (a+b \log \left (c x^n\right )\right )}{3 b n} \]

output
2/3*sec(a+b*ln(c*x^n))^(3/2)*sin(a+b*ln(c*x^n))/b/n+2/3*(cos(1/2*a+1/2*b*l 
n(c*x^n))^2)^(1/2)/cos(1/2*a+1/2*b*ln(c*x^n))*EllipticF(sin(1/2*a+1/2*b*ln 
(c*x^n)),2^(1/2))*cos(a+b*ln(c*x^n))^(1/2)*sec(a+b*ln(c*x^n))^(1/2)/b/n
 
3.3.71.2 Mathematica [A] (verified)

Time = 0.16 (sec) , antiderivative size = 69, normalized size of antiderivative = 0.74 \[ \int \frac {\sec ^{\frac {5}{2}}\left (a+b \log \left (c x^n\right )\right )}{x} \, dx=\frac {2 \sec ^{\frac {3}{2}}\left (a+b \log \left (c x^n\right )\right ) \left (\cos ^{\frac {3}{2}}\left (a+b \log \left (c x^n\right )\right ) \operatorname {EllipticF}\left (\frac {1}{2} \left (a+b \log \left (c x^n\right )\right ),2\right )+\sin \left (a+b \log \left (c x^n\right )\right )\right )}{3 b n} \]

input
Integrate[Sec[a + b*Log[c*x^n]]^(5/2)/x,x]
 
output
(2*Sec[a + b*Log[c*x^n]]^(3/2)*(Cos[a + b*Log[c*x^n]]^(3/2)*EllipticF[(a + 
 b*Log[c*x^n])/2, 2] + Sin[a + b*Log[c*x^n]]))/(3*b*n)
 
3.3.71.3 Rubi [A] (verified)

Time = 0.36 (sec) , antiderivative size = 91, normalized size of antiderivative = 0.98, number of steps used = 8, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.368, Rules used = {3039, 3042, 4255, 3042, 4258, 3042, 3120}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\sec ^{\frac {5}{2}}\left (a+b \log \left (c x^n\right )\right )}{x} \, dx\)

\(\Big \downarrow \) 3039

\(\displaystyle \frac {\int \sec ^{\frac {5}{2}}\left (a+b \log \left (c x^n\right )\right )d\log \left (c x^n\right )}{n}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\int \csc \left (a+b \log \left (c x^n\right )+\frac {\pi }{2}\right )^{5/2}d\log \left (c x^n\right )}{n}\)

\(\Big \downarrow \) 4255

\(\displaystyle \frac {\frac {1}{3} \int \sqrt {\sec \left (a+b \log \left (c x^n\right )\right )}d\log \left (c x^n\right )+\frac {2 \sin \left (a+b \log \left (c x^n\right )\right ) \sec ^{\frac {3}{2}}\left (a+b \log \left (c x^n\right )\right )}{3 b}}{n}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\frac {1}{3} \int \sqrt {\csc \left (a+b \log \left (c x^n\right )+\frac {\pi }{2}\right )}d\log \left (c x^n\right )+\frac {2 \sin \left (a+b \log \left (c x^n\right )\right ) \sec ^{\frac {3}{2}}\left (a+b \log \left (c x^n\right )\right )}{3 b}}{n}\)

\(\Big \downarrow \) 4258

\(\displaystyle \frac {\frac {1}{3} \sqrt {\sec \left (a+b \log \left (c x^n\right )\right )} \sqrt {\cos \left (a+b \log \left (c x^n\right )\right )} \int \frac {1}{\sqrt {\cos \left (a+b \log \left (c x^n\right )\right )}}d\log \left (c x^n\right )+\frac {2 \sin \left (a+b \log \left (c x^n\right )\right ) \sec ^{\frac {3}{2}}\left (a+b \log \left (c x^n\right )\right )}{3 b}}{n}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\frac {1}{3} \sqrt {\sec \left (a+b \log \left (c x^n\right )\right )} \sqrt {\cos \left (a+b \log \left (c x^n\right )\right )} \int \frac {1}{\sqrt {\sin \left (a+b \log \left (c x^n\right )+\frac {\pi }{2}\right )}}d\log \left (c x^n\right )+\frac {2 \sin \left (a+b \log \left (c x^n\right )\right ) \sec ^{\frac {3}{2}}\left (a+b \log \left (c x^n\right )\right )}{3 b}}{n}\)

\(\Big \downarrow \) 3120

\(\displaystyle \frac {\frac {2 \sin \left (a+b \log \left (c x^n\right )\right ) \sec ^{\frac {3}{2}}\left (a+b \log \left (c x^n\right )\right )}{3 b}+\frac {2 \sqrt {\sec \left (a+b \log \left (c x^n\right )\right )} \sqrt {\cos \left (a+b \log \left (c x^n\right )\right )} \operatorname {EllipticF}\left (\frac {1}{2} \left (a+b \log \left (c x^n\right )\right ),2\right )}{3 b}}{n}\)

input
Int[Sec[a + b*Log[c*x^n]]^(5/2)/x,x]
 
output
((2*Sqrt[Cos[a + b*Log[c*x^n]]]*EllipticF[(a + b*Log[c*x^n])/2, 2]*Sqrt[Se 
c[a + b*Log[c*x^n]]])/(3*b) + (2*Sec[a + b*Log[c*x^n]]^(3/2)*Sin[a + b*Log 
[c*x^n]])/(3*b))/n
 

3.3.71.3.1 Defintions of rubi rules used

rule 3039
Int[u_, x_Symbol] :> With[{lst = FunctionOfLog[Cancel[x*u], x]}, Simp[1/lst 
[[3]]   Subst[Int[lst[[1]], x], x, Log[lst[[2]]]], x] /;  !FalseQ[lst]] /; 
NonsumQ[u]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3120
Int[1/Sqrt[sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2/d)*EllipticF[(1/2 
)*(c - Pi/2 + d*x), 2], x] /; FreeQ[{c, d}, x]
 

rule 4255
Int[(csc[(c_.) + (d_.)*(x_)]*(b_.))^(n_), x_Symbol] :> Simp[(-b)*Cos[c + d* 
x]*((b*Csc[c + d*x])^(n - 1)/(d*(n - 1))), x] + Simp[b^2*((n - 2)/(n - 1)) 
  Int[(b*Csc[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1] 
&& IntegerQ[2*n]
 

rule 4258
Int[(csc[(c_.) + (d_.)*(x_)]*(b_.))^(n_), x_Symbol] :> Simp[(b*Csc[c + d*x] 
)^n*Sin[c + d*x]^n   Int[1/Sin[c + d*x]^n, x], x] /; FreeQ[{b, c, d}, x] && 
 EqQ[n^2, 1/4]
 
3.3.71.4 Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(290\) vs. \(2(119)=238\).

Time = 20.80 (sec) , antiderivative size = 291, normalized size of antiderivative = 3.13

method result size
derivativedivides \(-\frac {2 \left (-2 \sqrt {\frac {1}{2}-\frac {\cos \left (a +2 b \ln \left (\sqrt {c \,x^{n}}\right )\right )}{2}}\, \sqrt {-1+2 {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}}\, \operatorname {EllipticF}\left (\cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right ), \sqrt {2}\right ) {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}-2 {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2} \cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )+\sqrt {\frac {1}{2}-\frac {\cos \left (a +2 b \ln \left (\sqrt {c \,x^{n}}\right )\right )}{2}}\, \sqrt {-1+2 {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}}\, \operatorname {EllipticF}\left (\cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right ), \sqrt {2}\right )\right ) \sqrt {\left (2 {\cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}-1\right ) {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}}}{3 n \sqrt {-2 {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{4}+{\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}}\, {\left (2 {\cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}-1\right )}^{\frac {3}{2}} \sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right ) b}\) \(291\)
default \(-\frac {2 \left (-2 \sqrt {\frac {1}{2}-\frac {\cos \left (a +2 b \ln \left (\sqrt {c \,x^{n}}\right )\right )}{2}}\, \sqrt {-1+2 {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}}\, \operatorname {EllipticF}\left (\cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right ), \sqrt {2}\right ) {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}-2 {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2} \cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )+\sqrt {\frac {1}{2}-\frac {\cos \left (a +2 b \ln \left (\sqrt {c \,x^{n}}\right )\right )}{2}}\, \sqrt {-1+2 {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}}\, \operatorname {EllipticF}\left (\cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right ), \sqrt {2}\right )\right ) \sqrt {\left (2 {\cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}-1\right ) {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}}}{3 n \sqrt {-2 {\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{4}+{\sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}}\, {\left (2 {\cos \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right )}^{2}-1\right )}^{\frac {3}{2}} \sin \left (\frac {a}{2}+\frac {b \ln \left (c \,x^{n}\right )}{2}\right ) b}\) \(291\)

input
int(sec(a+b*ln(c*x^n))^(5/2)/x,x,method=_RETURNVERBOSE)
 
output
-2/3/n*(-2*(sin(1/2*a+1/2*b*ln(c*x^n))^2)^(1/2)*(-1+2*sin(1/2*a+1/2*b*ln(c 
*x^n))^2)^(1/2)*EllipticF(cos(1/2*a+1/2*b*ln(c*x^n)),2^(1/2))*sin(1/2*a+1/ 
2*b*ln(c*x^n))^2-2*sin(1/2*a+1/2*b*ln(c*x^n))^2*cos(1/2*a+1/2*b*ln(c*x^n)) 
+(sin(1/2*a+1/2*b*ln(c*x^n))^2)^(1/2)*(-1+2*sin(1/2*a+1/2*b*ln(c*x^n))^2)^ 
(1/2)*EllipticF(cos(1/2*a+1/2*b*ln(c*x^n)),2^(1/2)))*((2*cos(1/2*a+1/2*b*l 
n(c*x^n))^2-1)*sin(1/2*a+1/2*b*ln(c*x^n))^2)^(1/2)/(-2*sin(1/2*a+1/2*b*ln( 
c*x^n))^4+sin(1/2*a+1/2*b*ln(c*x^n))^2)^(1/2)/(2*cos(1/2*a+1/2*b*ln(c*x^n) 
)^2-1)^(3/2)/sin(1/2*a+1/2*b*ln(c*x^n))/b
 
3.3.71.5 Fricas [C] (verification not implemented)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.08 (sec) , antiderivative size = 145, normalized size of antiderivative = 1.56 \[ \int \frac {\sec ^{\frac {5}{2}}\left (a+b \log \left (c x^n\right )\right )}{x} \, dx=\frac {-i \, \sqrt {2} \cos \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right ) {\rm weierstrassPInverse}\left (-4, 0, \cos \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right ) + i \, \sin \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right )\right ) + i \, \sqrt {2} \cos \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right ) {\rm weierstrassPInverse}\left (-4, 0, \cos \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right ) - i \, \sin \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right )\right ) + \frac {2 \, \sin \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right )}{\sqrt {\cos \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right )}}}{3 \, b n \cos \left (b n \log \left (x\right ) + b \log \left (c\right ) + a\right )} \]

input
integrate(sec(a+b*log(c*x^n))^(5/2)/x,x, algorithm="fricas")
 
output
1/3*(-I*sqrt(2)*cos(b*n*log(x) + b*log(c) + a)*weierstrassPInverse(-4, 0, 
cos(b*n*log(x) + b*log(c) + a) + I*sin(b*n*log(x) + b*log(c) + a)) + I*sqr 
t(2)*cos(b*n*log(x) + b*log(c) + a)*weierstrassPInverse(-4, 0, cos(b*n*log 
(x) + b*log(c) + a) - I*sin(b*n*log(x) + b*log(c) + a)) + 2*sin(b*n*log(x) 
 + b*log(c) + a)/sqrt(cos(b*n*log(x) + b*log(c) + a)))/(b*n*cos(b*n*log(x) 
 + b*log(c) + a))
 
3.3.71.6 Sympy [F(-1)]

Timed out. \[ \int \frac {\sec ^{\frac {5}{2}}\left (a+b \log \left (c x^n\right )\right )}{x} \, dx=\text {Timed out} \]

input
integrate(sec(a+b*ln(c*x**n))**(5/2)/x,x)
 
output
Timed out
 
3.3.71.7 Maxima [F]

\[ \int \frac {\sec ^{\frac {5}{2}}\left (a+b \log \left (c x^n\right )\right )}{x} \, dx=\int { \frac {\sec \left (b \log \left (c x^{n}\right ) + a\right )^{\frac {5}{2}}}{x} \,d x } \]

input
integrate(sec(a+b*log(c*x^n))^(5/2)/x,x, algorithm="maxima")
 
output
integrate(sec(b*log(c*x^n) + a)^(5/2)/x, x)
 
3.3.71.8 Giac [F(-1)]

Timed out. \[ \int \frac {\sec ^{\frac {5}{2}}\left (a+b \log \left (c x^n\right )\right )}{x} \, dx=\text {Timed out} \]

input
integrate(sec(a+b*log(c*x^n))^(5/2)/x,x, algorithm="giac")
 
output
Timed out
 
3.3.71.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\sec ^{\frac {5}{2}}\left (a+b \log \left (c x^n\right )\right )}{x} \, dx=\int \frac {{\left (\frac {1}{\cos \left (a+b\,\ln \left (c\,x^n\right )\right )}\right )}^{5/2}}{x} \,d x \]

input
int((1/cos(a + b*log(c*x^n)))^(5/2)/x,x)
 
output
int((1/cos(a + b*log(c*x^n)))^(5/2)/x, x)